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Distortions in a family of conjugated polymers are stud- 
ied within two complementary approaches, i.e. within a many- 
body Valence Bond (VB) approach using a transfer matrix 
technique to treat the Heisenberg model of the systems, and 
also in terms of the tight-binding band-theoretic model with 
interactions limited to nearest neighbors. The computations 
indicate that both methods predict the presence or absence 
of the same distortions in most of the polymers studied. 
PACS 71.10 -fx, 71.27 -fa, 63.20.Kr, 31.20.Tz 



I. INTRODUCTION 

The recent discovery of the first inorganic spin-Peierls 
material, CuGeOs, ||l| has engendered a renewed inter- 
est in Spin-Peierls systems, i.e. systems which present a 
structural distortion below the spin-Peierls temperature 
due to residual magnetoelastic couplings stabilizing the 
ground-state, in analogy to Peierls distortion |^ asso- 
ciated to an electron-soft phonon instability opening a 
band gap at the Fermi level. Recent experiments [D sug- 
gest that this is not an isolated case, and the pronounced 
decrease of susceptibility observed |^ in Q;-NaV205 is 
also due to a spin-Peierls transition. 

The spin-Peierls transition was first observed in pre- 
dominantly organic compounds as TTFCuBDT ig, TTF- 
TCNQ ||, (TMTSF)2PF6 §, or TTF-AuBDT [|. 
Theoretically, it has been studied (see for instance 
|,|0|jll]jl|,|^,|l|l and Refs. therein) as a geometrical 
symmetry breaking for the lowest eigenstate of a Heisen- 
berg Hamiltonian. Peierls and spin-Peierls phenomena 
are still a subject of discussion for many other poly- 
mers, since if a deviation occurs that lowers the chain's 
symmetry, then different symmetry-equivalent distorted 
ground-states may arise which correspond to different 
thermodynamic phases and, at sufficiently low temper- 
ature, the possibility of solitonic excitations and/or con- 
duction could arise [ p^ . 

Furthermore, it has been argued [ p^ that under sim- 
ilar structural circumstances a Peierls distortion is pre- 
dicted for the simple Hiickel tight-binding model of tt- 
network strips if and only if a spin-Peierls distortion 
is also predicted from Valence Bond (VB) theory (or 
the formally equivalent s = 1/2 Heisenberg model) at 
the simple resonance theoretic level. At this level of 
approximation the VB wave-functions are restricted to 
equally weighted superpositions of special covalent VB 
singlet states, i.e. of Kekule structures |17|, where ev- 
ery TT-electron is coupled to a singlet state with one of 
their nearest-neighbors. These Kekule structures may be 
partitioned into long-range ordered spin-pairing phases, 
the lowest-lying phase corresponding to the highest count 
of Kekule structures contributing to it. Within this 
approach, a spin-Peierls distortion is predicted if there 
are two maximum-cardinality-degenerate Kekule phases 
(see ll8| and references therein). Then, this correspon- 
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dence between Peierls and spin-Peierls instabilities im- 
plies that a zero- width band gap for a 7r-network polymer 
is predicted if and only if there are two such cardinality- 
degenerate Kekule phases. The question then arises as 
to whether this correspondence is maintained when going 
beyond the resonance theoretic approximation. 

For instance, the dimerization in polyacetylene has 
traditionally been interpreted in terms of band the- 
ory [^,^ as a Peierls distortion. Recently, however, 
this dimerization has also been successfully explained 
pl| with a Heisenberg-like Hamiltonian model |2^] as 
a spin-Peierls distortion, using both cluster-expanded 
wave- functions and perturbation theory. This cluster- 
expanded many-body treatment of distortions has also 
been applied to the polyacene polymer j2^, which ear- 
lier has been extensively studied from the independent 
particle point of view, since it exhibits an accidental zero- 
width band gap at a simple tight-binding level (see [ p3|j24| 
and references therein), and a new quasi-degeneracy has 
been predicted. 

Thence, the comparison between the independent- 
particle and the many-body VB treatments for degener- 
acy and symmetry-breaking in polymers deserves further 
analysis. It is our purpose here to investigate the ground- 
state symmetries and degeneracies for several conjugated 
polymers using both a simple many-body VB framework 
and a simple tight-binding model. The rationale for these 
simplest models (with just nearest-neighbor interactions) 
is that they reveal distortive responses which qualitative 
dominate over the otherwise harmonic responses (e.g., 
associated with the a electrons. That is, these sim- 
plest models should reveal dominant qualitative features, 
which should persist independently of parameterization. 

The polymers we focus our attention on are: poly- 
aceacene (PAA), poly(benz[m,n])anthracene (PBA), and 
polyperylene (PPR) (Fig. ^ ) . All these systems exhibit a 
zero-width band gap at the simplest tight-binding level. 
So far, very few experimental results are available. Only 
PBA and PPR |^ have already been synthesized. 
Theoretically PPR has been treated from the indepen- 
dent particle point of view |2^j2^, and also using the 
valence effective hamiltonian technique ||2^ , while as far 
as we know PBA has not been previously treated. PAA 
has been discussed in the literature, mostly from an in- 
dependent particle point of view [p8| , ^()pl|j3^ , |3^ and less 
frequently from a resonance-theoretic approach |]l5| , |l6| , 
though it has not been synthesized yet. PAA can be 
seen, together with polyacetylene and polyacene, as the 
first members of a family of poly-trans-polyacetylenes, 
graphite being the final member of the family. All these 
can be thought of as special cases of ladder materials ||3^] 
as already had been pointed out in Refs. |l5| , |3^ , ^ . 

Within the many-body VB framework, we will con- 
sider the antiferromagnetically-signed spin-^ Heisen- 
berg model (for more general derivations of this model 
than those based on degenerate perturbation expansions 



see for instance |2^,^ and references therein). Ade- 
quate many-body wave-function ansdtze provide varia- 
tional upper bounds to the ground-state energy. Two 
different kinds of variational localized-site cluster ex- 
panded ansdtze have been considered: first a Resonat- 
ing VB (RVB) ansatz, where the trial wave-function is 
a weighted superposition over all singlets constructed as 
products of singlet pairs each involving two (not nec- 
essarily nearest-neighbor) sites at a time; and second a 
Neel-state-based ansatz, where a Neel state is the zeroth 
order wave-function from which the trial wave-function 
is generated. We evaluate the matrix elements for each 
ansatz with a transfer-matrix technique introduced pre- 
viously [|l],|3|j3^,|^,||,|§. For the tight-binding band 
theory calculations we consider the so-called translation- 
ally adapted Hiickel model limited to nearest neighbors. 

This paper is organized as follows: in section || the de- 
scription of the polymers, their symmetries and relevant 
distortions are given. In section III we introduce briefly 
the translationally adapted Hiickel model. In section IV 
a description of the VB method is given in terms of the 
Heisenberg Hamiltonian and the trial wave-functions are 
presented. Also, the technique to compute the physical 
magnitudes based on a transfer matrix is introduced and 
applied to obtain the ground-state energy of the systems. 
Results are presented and discussed in section |v|. And 



finally our conclusions can be found in section VI 



II. DESCRIPTION OF THE POLYMERS AND 
THEIR SYMMETRIES 

The systems studied are polymeric strips of finite width 
and infinite length (L oo) (see Fig. |l|). They are con- 
structed with fused benzene rings, and can be seen as 
cut from the 2-dimensional graphite or honeycomb lat- 
tice. Each site of the lattice is taken to represent an 
sp^-hybridized carbon atom with one 7r-orbital perpen- 
dicular to the plane of the lattice and with one 7r-electron 
per site. These strips are presumed to be translation- 
ally symmetric along L, with periodic boundary condi- 
tions; so that the strips may be divided into unit cells 
or eventually reduced unit cells, when the space group 
of the strip contains operations involving glide reflec- 
tions. The space group of the strips include, along with 
the primitive translation, rotations C„, reflections a, and 
combination of rotations and reflections (improper rota- 
tions), coordinate inversion i, and screw-rotations and 
glide-reflections, Cs, i.e. a combination of an improper 
two-fold rotation or reflection with a nonprimitive trans- 
lation of half a unit cell which by themselves do not leave 
the lattice invariant (see for instance Fig. ^ and Table |). 

Of special interest are minimal subsets of symmetry 
operations, whose removal lead to: (i) a band gap open- 
ing at the Fermi level, when analyzed from the band- 
theoretic point of view; (ii) the lifting of the degeneracy 
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of Kekule phases, if seen from the resonance-theoretic 
treatment. If a zero gap occurs at fc = tt (as is frequently 
the case for benzenoid polymers) then such a minimal 
subset will be so as to no more than double the size of a 
unit cell. 

When a symmetry is broken, there is a distortion pa- 
rameter A; associated to the stretching or shortening of 
the bond I, with bonds numbered as in Figs. |, I and 
^. Two symmetry elements are chosen to label the inter- 
esting distortions for every polymer as shown in Tables 
||, m and IV, where appropriate constraints on the A;, 
imposed by the different symmetry breakings, are also 
shown. The distortions are classified as to symmetric 
or antisymmetric (-1) with respect to these two se- 
lected symmetry elements. 

The PAA polymer is formed by benzene rings sharing 
four consecutive edges with neighboring rings as shown in 
Fig. 0(a). It can also be seen as a trimer of non-dimerized 
parallel all-trans polyacetylene chains. The 6-site unit 
cell can be broken into two 3-site reduced unit cells, de- 
fined as the region between dashed lines in Fig. |l|(a). In 
the band picture, there is a half-filled band and, conse- 
quently, a zero-width band gap is predicted, regardless 
of distortions which preserve the glide-reflection symme- 
try. In the simplest VB picture, i.e. resonance theory, 
there are two maximum-cardinality degenerate Kekule 
phases. For instance, defining M as the number of "dou- 
ble bonds" crossed by an oblique line (see Fig. there 
are two Kekule phases M— even equivalent to two M— 
odd which do not mix because of the cyclic boundary 
conditions of the strip and they are degenerate since they 
each contain essentially a single Kekule structure. A dis- 
tortion that could open the band gap at the Fermi level 
and lift the degeneracy of the Kekule phases requires the 
destruction of the glide-refiection symmetry. The distor- 
tions to be considered are then those which are antisym- 
metric with respect to interchange of the two reduced 
unit cells in a new unit cell. 

PBA is formed by a polyacene strip where added ben- 
zene rings have been, top and bottom, alternatively fused 
on (see Fig. 0(b)). A reduced unit cell can be de- 
fined for this system between the dashed lines in Fig. 
|l|(b), with two 7-sites reduced unit cells per unit cell. 
It is a half-filled band system and, like PAA, a zero- 
width band gap is predicted. Resonance theory, following 
Ref. fl^ , predicts two maximum-cardinality degenerate 
Kekule phases. As in the PAA polymer, the interest- 
ing distortions that could open the band gap and lift the 
degeneracy are those that are antisymmetric under op- 
erations which interchange the two types of reduced unit 
cells. 

The PPR polymer is formed by fused benzene rings 
as drawn in Fig. 0(c). The miit cell containing 10 sites 
is defined between the dashed lines in the graph, and 
there is no smaller reduced unit cell for this system. The 
space group is generated by the point group D2h and 



the translation operations along the strip (see Table | 
and Fig. ||(c)). Differently from the rest of the polymers 
here, there is no glide-reflection symmetry operation for 
PPR. Furthermore, it does not have an odd number of 
TT-electrons per reduced unit cell so that it does not cor- 
respond to a half-filled band system. Nevertheless, there 
is an accidental degeneracy at the Hiickel level of ap- 
proximation, so that it has a zero-width band gap any- 
way (see next section). Correspondingly, resonance the- 
ory predicts two maximum-cardinality degenerate Kekule 
phases p6|. A totally symmetric distortion will also be 



considered for this system (see Table IV) 



III. TRANSLATIONALLY ADAPTED HUCKEL 
MODEL 



The Hiickel model is the simplest tight-binding model: 

Huuck = y ^ I3ni,mj (c^io-Cmjo- + C.-nija'^niij) (1) 



E 

<ni,mj>,cr 



•^nio- (cnio-) are the creation (annihilation) electron op- 
erators on site i of unit cell n with spin a and Pni,mj 
is the "Hiickel resonance integral" (or hopping integral) 
between sites i and j in unit cells n and m, respectively. 
< ni,mj > indicates that the sum is restricted to near- 
est neighbors. Considering the translational invariance 
symmetry of the system, we can define translationally 
symmetry adapted states 
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^ n— 1 



(2) 



The matrix elements of the Hamiltonian between these 
new states are 



{j;k\H\i;k')^Skk' 



^ ^ ^—ik{n—m) 
<ni,mj'> 



f3n 



(3) 



Diagonalizing the hamiltonian matrix elements, the en- 
ergy bands e{k) are finally obtained. 

Symmetry breaking can be considered taking (3ni,mj as 



f3n 



/3(1 + A, 



(4) 



where A„i^mj (| ^ni,mj 1^ 1) is the distortion parameter, 
as introduced in the previous section, that measures the 
strength of the distortion between sites ni and mj. 

IV. VALENCE BOND METHOD 

Within the VB picture we attempt here to go beyond 
resonance theory when solving the Heisenberg Hamilto- 
nian: 
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Hhcis — / ^ Jni.rajSniSraj (5) 

Jni,mj is the "exchange integral" between nearest- 
neighbor sites ni and mj and Sni denotes the spin oper- 
ator for one electron on site ni. 

Jni.rnj — ^ni^rnj^ (6) 

with Ani^mj being the distortion parameter associated 
to the bond between sites ni and mj when there is a 
symmetry breaking. 

While solving the Hiickel model is an easy task, solv- 
ing the Heisenberg hamiltonian is in general a non-trivial 
problem. In order to obtain, along with the appropri- 
ate approximate wave- functions, good variational upper 
bounds to the ground-state energy of this model, E{A), 
for the polymer systems, we consider two different types 
of cluster-expanded ansdtze that depend on variational 
parameters, each of which describes the local features of 
the system. Since our polymers are bipartite systems 
with total spin zero, we have considered a Neel-state- 
based ansatz and a RVB ansatz. These ansdtze were 
introduced in Ref. |23) and we shall make here a brief 
description of them. Related ansdtze have also been suc- 
cessfully considered by other authors |^,^ when solv- 
ing the s = 1/2 Heisenberg hamiltonian for the square 
lattice. 

A. Neel-state-based ansatz (NSBA) 

The cluster expanded wave-function ansatz in this sec- 
tion is based upon the Neel state as a zeroth-order wave- 
function, 

ieA j&B 

I = n n '^(^■) 

i j 

where A and B denote the two sets of sites such that each 
member of one set is a nearest neighbor solely to (some) 
sites of the other set, and a{i) (/?(«)) indicate that the 
spin of the electron on site i is +1/2 (—1/2). A lowering 
of the energy, with respect to that of the Neel state, oc- 
curs for an ansatz defined within a subspace spanned by 
I <i>jv) and the states obtained when applying to | $jv) 
the XY terms, S^-S^j, of the Heisenberg operator, an 
arbitrary number of times in an "unlinked" way. These 
additional states which are to be mixed with the Neel 
state can be generated in terms of the nearest-neighbor 
pair excitation operator: 

GA <ni,mj> 
= 5Z ^ ^ni,m]S,,^^S^j (8) 
ni mj 

where the Xni,mj are scalars to be optimized and S^^^ and 
•^mi spm raising and lowering operators on site ni 



St,^Sl,±\Sl^ (9) 

From that, the Neel-state-based wave- function ansatz 
(NSBA) will be a cluster-expansion in terms of P excita- 
tions acting on the Neel state 

I *Ar) = f/e^ I (10) 

where U indicates that only unlinked terms are to be re- 
tained from the Taylor series expansion. Namely, | ^ n) 
is a wave-function where the Neel state is mixed with 
states that differ from it by an arbitrary number of cou- 
ples of disjoint pairs of neighboring spins that have been 
flipped, each state in the superposition being weighted 
by the product of the variational parameters associated 
to the flips in that state. 

B. Resonating Valence Bond ansatze 

In this approach we start with a one-bond-range 
RVB (IBR-RVB), that plays the fundamental zeroth- 
order role for the more elaborated three-bond-range RVB 
(3BR-RVB) ansatz following. 

1. One- Bond- Range RVB ansatz 

A IBR-RVB I ^i) is a weighted superposition of 
Kekule states, i.e. nearest neighbor VB states, where ev- 
ery site ni is spin paired to one of its neighbors mj. It 
can be written as 

GA <ni,mj> 

I *1> = C^O n E - ^n^Stai) I n) (H) 

ni mj 

I is the identity operator, Uq indicates that the terms to 
be retained are those where each site appears once and 
only once, and the weighting factor of a Kekule state 
in ^E"! is a product of variational parameters Xni,mj as- 
sociated to the singlet pairs ni,mj in the Kekule state 
considered. 

2. Three- Bond- Range RVB ansatz 

The 3BR-RVB is a weighted superposition of all the 
VB structures within a phase with each spin-pairing be- 
tween A- and B-sublattice sites separated by no more 
than 3 bonds. In the usual form for cluster expanded 
wave- functions, it may be viewed as generated from the 
IBR-RVB as follows: The XY terms, S^^S^^, of the 
Heisenberg Hamiltonian acting on the IBR-RVB wave- 
function of Eqn. (|ri|) yield "long-bonded" states with 
pairings among three-bond distant neighbors, along with 
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"neighbor-bonded" states already in V&i. These "long- 
bonded" states can be directly generated by the "recou- 
pling" of two neighboring bond-singlets in ^'i (see Ref. 
psf). We may denote by qef the operator related to such 
a recoupling between two bond-singlets e and /. From 
I ^i) we may build the 3BR-RVB allowing an arbitrary 
number of recouplings of two simply neighboring bond- 
singlets, i.e. unlinked pairs with one and only one site in a 
pair being a nearest neighbor to a site in the other pair. 
Then, the overall 3BR-RVB excitation operator above 
the IBR-RVB wave-function might be viewed to be 



J2 

<e,f> 



(12) 



with the Xef being variational parameters, and where 
< e, / > indicates that the sum is restricted to sim- 
ply neighboring bond-singlets. The corresponding ansatz 
would then be 



[/e'3 I ^fi) 



(13) 



where again U indicates that only unlinked terms are to 
be retained. That is, in the Taylor series expansion of 
e*^ one retains only products of q^f such that no pair 
index (e or /) shares any vertices with another pair in- 
dex in the product. And Q and are to be optimized 
simultaneously. 



Tts = {et\T \ es) 



(15) 



denote a weighted sum over the various ways a local state 
Cs may succeed a local state ej. The weight of every 
contribution is obtained by considering the variational 
parameters associated to the way et evolves to Cs, and, 
eventually, additional factors coming from Pauling's su- 
perposition rules The overlap is then evaluated in 
terms of the T matrix: 



(* I *) 



trT^ 



(16) 



For L ^ oo, the largest eigenvalue A of T dominates, 
and the overlap reduces to 



(* I 



A^ 



(17) 



The hamiltonian expectation value over | can be 
obtained in a similar way introducing a "connection" ma- 
trix, C, defined according to 

per cell 



<ini,mj> 



(18) 



where c measures the range of the interaction within the 
ansatz. In our case c = 2. And, the matrix element 



C. Expectation- value calculations by the 
Transfer-Matrix Technique 



The ground-state energy 



(14) 



is computed as a function of variational parameters for 
each of the above introduced ansdtze assuming transla- 
tional symmetry and cyclic boundary conditions along L. 
The way our ansdtze are chosen allows us to deal with 
the systems locally, so that one can define a transfer ma- 
trix j23j that describes the local features and reduces the 
computation of Eqn. (|^ to products of "small" matrices 
pl| , p3 36 3^,^,^,^. Let's suppose there are imaginary 
vertical lines cutting the strip on translationally equiva- 
lent positions (including improper translations). We can 
define the a^sate-dependent "local states" according to 
every possible local spin-pairing/spin-flip pattern around 
a given position determined by one of the imaginary ver- 
tical lines, and ultimately use this to compute | ^). 
Thence these local states contain the contributions from 
both the bra and the ket. From the assumed transla- 
tional symmetry, local states in every position are to be 
the same. Now, labelling these local states by et, t rang- 
ing over the whole set of local states, we let the transfer- 
matrix element 



C'ts = (et I C I e,) 



(19) 



is a weighted sum over the various ways a local state Cs 
may succeed a local state et after c transfer-matrix-steps 
when the Hamiltonian operators per unit cell are present. 
In the long-length limit Eqn. (|lj) reduces to: 



E = 



1 (A,/|C| A,r) 
A2 (A,;|A,r) 



(20) 



where {A, I \ and | A,r) are left and right eigenvectors 
corresponding to the maximum eigenvalue A of T. This 
expression is a function of the variational parameters as- 
sociated to 5* and an upper bound to the exact ground- 
state energy is obtained. Implementation of a suitable 
numerical optimization yields a best upper bound. The 
energy expression ( |20| ) can be readily generalized when 
considering possible distortions. The connection matrix 
per unit cell can be understood as a sum of matrices 
Cni,mj, each one concerning two body interactions be- 
tween neighboring sites ni and mj, weighted by the factor 
ni.mj that modifies its interaction strength. Then 



1 



Gn <ni,mj> 
i mj 



(21) 



and the energy expectation evaluation is reduced to some 
"simple" matrix manipulations. 
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V. RESULTS AND DISCUSSION 

Computations based on band theory at a Hiickel tight- 
binding level of approximation (see section [II ), and 
within VB theory with the cluster expanded IBR-RVB 
ansatz and NSBA (see section IV) were carried out for 
all polymer systems here described. For the PAA also 
the 3BR-RVB has been used. In this case the IBR-RVB 
ansatz contains only one Kekule structure, so it is es- 
pecially appropriate to use the 3BR-RVB wave-function 
and go beyond a single Kekule-structure approximation. 
This circumstance differs with the rest of the polymers, 
where the number of Kekule structures in the correspond- 
ing IBR-RVB is big. The different VB upper bounds 
to the energy of the undistorted polymers are presented 
in Table ^ together with that for polyacetylene. The 
lowest upper bound to the ground-state energy for the 
undistorted system is given by the NSBA. 



1. PAA 

The highest occupied Hiickel tight-binding band and 
the lowest unoccupied band cross at fc = tt. Taking 

in the Hiickel reso- 



into account the perturbation 



nance integral, P„ 



(3(1 



), only the totally- 



antisymmetric distortion with respect to Cs and la- 
belled as C3 (see Table ||) opens a gap at k = tt. But 
the leading term of the energy lowering AE versus A is 
~ A^, as it is the positive phonon energy contribution 
to be added. So band theory at this low level of ap- 
proximation predicts neither a presence nor absence of 
a C3 distortion for this system, the result depending on 
the final balance between these two contributions to the 
energy. Nevertheless, if interactions with more distant 
TT-centers are included, although small, linear terms in A 
are argued to arise |16| and then the distortion is favored. 

Still within band theory, this system has also been 
studied by other authors at different levels of approxi- 
mation. Kertesz and Tanaka |^ suggest a totally 
antisymmetric distortion, though leading to a quadratic 
small gap that could be suppressed by interchain interac- 
tions. The tight-binding SCF-MO method at the level of 
CNDO/2 (complete neglect of differential overlap) cal- 
culations suggests that the Peierls distortion does not 
take place so one can expect metallic behavior ||30|| , while 
Bozovic ||2^ combining tight-binding calculations with 
group-theoretical arguments predicts distortions of type 
B (see Table ||) as favored. Therefore, within band the- 
ory, predictions about the opening or not of a band gap 
at the Fermi level, or the distortion driving it, depend 
crucially on the level of approximation. 

Let us consider now the many-body VB method. The 
ground-state energy has been obtained using the NSBA 
and both the IBR- and the 3BR-RVB ansdtze of sec- 
tion as a function of A for the different distortions 



A, i?i , i?2 , Ci , C2 , C3 (see Table ||) . Transfer and connec- 
tion matrices of dimensions 14 x 14 (for the NSBA) and 
60 X 60 (for the 3BR-RVB ansatz) were needed in order 
to carry out computations. The energy for the different 
distortions when the NSBA is used has been plotted as 
a function of A in Fig. 0(a), while results obtained with 
the 3BR-RVB ansatz are presented in Fig. 0(b). Plots 
from the IBR-RVB ansatz are not given, since they are 
qualitatively identical to those from the 3BR-RVB ones. 
Comparing NSBA and RVB ansdtze, it can be seen that 
the ordering of AE for the different distortions is the 
same in any case, the strongest lowering corresponding 
to the C3 distortion. 

Nevertheless, while the energy response to A, B and C 
distortions is linear for the RVB ansdtze, clearly predict- 
ing a C3 distorted ground-state, in the NSBA case they 
still go as ~ A^. Fitting the results in a parabolic curve, 
it is obtained that AE ~ — 1.923 A^. Again a distortion 
is not clearly predicted with our NSBA. A comparison 
of the coefficients coming from this term and those from 
the phonon energy should be made in order to decide 
whether this ansatz is able to predict or not to predict 
a C3 distortion. This ambiguity of prediction in some 
sense rationalizes earlier contradictory results: via nu- 
merical band theory by Yamabe ct al. | |30[ | , predicting an 
undistorted ground-state, and via band/group-theoretic 
considerations by Bozovic Q, predicting a B distortion. 

Although the RVB ground-state energy is higher than 
the NSBA, its predictions on ground-state instabilities 
are based upon the known giobal-singiet character of 
the ground state along with its local-singlct character, 
leading to asymptotically orthogonal and noninteracting 
phases responding essentially independently to distor- 
tions. Relaxation of this local-singlet character would im- 
ply the inclusion of pairings between distant sites, lead- 
ing to undesirable long-range correlations of the type in 
the Neel state. Then we expect a C3 distorted ground- 
state as predicted by RVB. Furthermore, NSBA at this 
lower level, with only two-site excitations, does not al- 
ways seem sensitive to instabilities as at higher order, 
such as for polyacetylene |^l|, then we expect that the 
distortion could be clearly predicted when going to a 
higher order NSBA too. Also it can be argued that in- 
clusion of slightly longer-range interactions (as between 
next-nearest neighbors) in the hamiltonian will increase 
the 'frustration' and the NSBA energy whereas the RVB 
expectations will change but little. Thus there is a ten- 
dency to invert the energy ordering of these states. Still 
another argument favoring RVB predictions is that the 
NSBA is not a pure singlet, as the ground state is known 
to be. Also the RVB type ansatz accords more closely to 
a classical organic chemical view of these polymers. 
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2. PBA 



VI. CONCLUSIONS 



The lowest occupied Hiickel tight-binding band and the 
highest unoccupied one cross ai k — n, so it is a zero- 
width band gap system. From aU the possible distortions 
considered in Table III, only the totally antisymmetric 



distortions Ci and C2 open a gap with an energy depen- 
dence linear in A. Therefore, band theory predicts that 
the system will distort. In the VB picture the possible 
distortions in PBA have been studied with the IBR-RVB 
ansatz. For this system we only carried out calculations 
with this ansatz for two reasons i) the IBR-RVB ansatz 
gives already a good upper bound to the ground-state 
energy because there is mixing of Kekule states and ii) 
the dimension of the transfer and connection matrices 
for the 3BR-RVB and the Neel-state-based ansdtze grow 
substantially with respect to the IBR-RVB one. In Fig. 
^ the energy of the IBR-RVB ansatz has been plotted 
as a function of A for the distortions Ai, A2, Ci and 
C2 classified in Table [II. The most favored distortions 
are the totally-antisymmetric ones Ci and C'2, in par- 
ticular Ci with a dependence ~ A. This result agrees 
with the predictions given from band theory, concluding 
that complementary approaches lead to the same kind of 
distortions for this system. 



3. PPR 

PPR is not a half-filled band system but the Hiickel 
model predicts an accidental zero-width band gap at fc = 
0. A, B, C and D distortions (see Table have been 
considered. The distortions Ci and C2 open a gap at fc = 
weakly, with an energy dependence AE ~ A^. But the 
totally symmetric distortion Di opens a band gap with 
an energy response linear in A. This result agrees with 



the predictions given by Bozovic |28| and Tanaka et al. 
prf . In Fig. ^(a) the Neel-state-based energy obtained, 
using 5x5 transfer and connection matrices, is plotted 
as a function of A for various possible distortions (see 



Table [V). Clearly the totally symmetric distortion, Di 
is favored with a linear energy dependence on A. Also in 
Fig. 11(b) the IBR-RVB energy is plotted for the various 
distortions as a function of A and results agree with the 
NSBA energy, namely that the Di distortion is the most 
favored one with a linear dependence in A. As in PBA, 
the IBR-RVB ansatz gives already a good upper bound 
due to the mixing of Kekule states. 

Band theory and the many-body VB method predict 
the same distortional behavior for this system, i.e. the 
system is unstable to a totally symmetric Di distortion. 
Some evidence exists for polyperylene synthesis [ p6| , but 
further experimental information on the structure (and 
properties) of this system is still needed. 



We have presented, both with the simple Hiickel 
tight-binding band theory and with a Heisenberg model 
hamiltonian (or, equivalently, the VB model), a study 
of the ground-state nature of a family of polymers: 
polyaceacene, poly(benz[m,n] anthracene), and polyper- 
ylene. We have focused our attention on correspon- 
dences between Peierls and spin-Peierls instabilities pre- 
dictions, when analyzed from these two complementary 
approaches. 

Upper bounds to the energy of the Heisenberg model 
in each case have been obtained with two alternative 
localized-site cluster-expanded wave-functions, i.e. RVB- 
type ansdtze and a Neel-based ansatz. We have shown 
that simple expressions of the physical magnitudes we 
were interested in, were easily obtained by using the 
transfer-matrix technique of ref. [p3| . 

From our results, it is concluded that the RVB wave- 
functions considered, which are restricted to IBR type for 
all the systems other than PAA, do not give our best up- 
per bound to the ground-state energy of the undistorted 
systems. Nevertheless, they are relevant for studying 
such phenomena as the spin-Peierls instability and el- 
ementary excitations as hole excitations or excitonic ex- 
citations as already pointed out [Q. Moreover the RVB 
ansdtze have a global-singlet character and a local-singlet 
character, precluding long-range order of the type of the 
Neel state, and generally improve relative to Neel-based 
ansdtze upon inclusion of higher-order (frustrative) terms 
in an elaborated Heisenberg model. 

The Neel-state-based ansatz gives a fairly good up- 
per bound to the ground-state energy for all the sys- 
tems considered. For the nearest-neighbor model con- 
sidered, this ansatz always yields lower energy than the 
RVB ones for undistorted systems. We have shown that 
with such a simple Neel-state-based wave-function, the 
corresponding energy is notably lower than the energy of 
the Neel state, while computations remain fairly simple. 
The Neel-state-based ansatz predicts for the polymers 
studied, the same distortions as the RVB description, 
except for the case of polyaceacene where this ansatz in 
our current simple considerations does not show whether 
the distortion is going to take place or not, although the 
strongest lowering of the energy also correspond a totally 
antisymmetric distortion. 

From the Heisenberg hamiltonian, or equivalently from 
the VB model, we have obtained that: 

1. PAA shows a totally antisymmetric distortion from 
the RVB, while NSBA is not conclusive, depending 
on the balance between electronic energy lowering 
and the phonon energy contribution. 

2. PBA shows a totally antisymmetric distortion. 



3. PPR is unstable to a totally symmetric distortion. 
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Within the band-theoretic picture, the Hiickel tight- 
binding model has been studied for all the same poly- 
mers. Results obtained for our 7r-network system are as 
follows: 

1. PAA could show a totally antisymmetric distortion 
at a simple Hiickel level, depending on the balance 
between electronic energy lowering and the phonon 
energy contribution. Other approximations already 
in the literature p8| , ^ yield contradictory results. 

2. PBA shows a totally antisymmetric distortion. 

3. PPR shows a totally symmetric distortion. 

Comparing band theory and the Heisenberg model 
results, it can be concluded that predictions of these 
two models based on opposite (or complementary) lim- 
its seem to lead to similar consequences under similar 
structural circumstances, i.e. both approaches predict the 
presence or absence of the same instability to symmetry 
for the polymers. It is to be noted that the band-theoretic 
results depend crucially on the level of approximation, as 
it is observed in the study of polyaceacene, where this 
picture at different levels of approximation gives rise to 
different predictions. On the other hand the Heisenberg 
model has proven to give predictions consistent from one 
level to another. Even in the case of PAA, where NSBA 
cannot make a clear prediction as it happens with band 
theory at its lower level, NSBA still shows the strongest 
lowering of the energy for the very same distortion sug- 
gested by RVB. Since the NSBA at this lower level, with 
only two-site excitations, does not always seems so sensi- 
tive to instabilities as at higher orders, such as it happens 
with polyacetylene |^. That is, the distortion some- 
times seems to only occur with a higher order NSBA, 
in agreement with RVB results. Therefore, it seems that 
the VB model, which includes correlation explicitly, gives 
a good description of these benzenoid systems, predict- 
ing spin-Peierls distortions whenever a Peierls distortion 
is also predicted. These results modify earlier sugges- 
tions (see ref. ' 44 1 ) that inclusion of correlation "a poste- 
riori", as a perturbation, diminishes the distortion. That 
is, we find any diminishment does not go to zero in the 
(strong correlation) Heisenberg-model limit, and indeed 
the RVB results indicate a stronger response to distor- 
tions (at least at the undistorted point on the potential 
energy hypersurface) . 

It has been shown that this treatment is computation- 
ally feasible specially for quasi-one dimensional systems 
where the transfer-matrix technique proves to be a pow- 
erful tool of computation. It is important to note that 
the results are developed in terms of quantities which 
remain finite as the strip-length goes to infinity. It is 
of some interest to compare the computational effort in- 
volved in the tight-binding approach versus that involved 
in our transfer-matrix cluster-expansion approach (for ei- 



ther RVB or Neel-state-based wavefunctions). The ma- 
trices H(fe) of eqn. I and T of eqn. |l^ arise in these 
respective approaches and are both finite independently 
of L — > oo. Both types are to be diagonalized, but there 
are some differences: 

1. Typically T increases in size much more rapidly 
with unit-cell "width" than does H(fc) (though 
these behaviours are reversed if the unit-cell 
"length" is considered instead). 

2. The total energy requires sampling many of the 
L — > oo H(fc) matrices (varying smoothly with 
wavevector k) whereas for given parameters there 
is but one T matrix to treat. 

3. The optimal total energy for the cluster expan- 
sions entails treating T matrices for numerous vari- 
ational parameter values whereas there is no much 
repetition with the H(fc). 

Notably if one goes beyond the tight-binding method to 
Hartree-Fock (or density-functional) approaches this last 
noted difference no longer occurs. Evidently the compu- 
tational effort via either SCF or our cluster expansion 
is roughly comparable (at least for linear polymers with 
modestly sized unit cells). 

The analysis carried out in this paper would require 
experimental testing. Though the synthesis of some of 
the systems considered like PAA seems quite difficult to 
achieve, there are hopes in this direction. Finally, some 
aspects of this treatment are not restricted only to the 
model Hamiltonian and the ground-state ansdtze pre- 
sented but can be applied to any system with effective 
short-range interactions if described by a localized-site 
cluster expanded ground-state wave- function. 
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TABLE I. Number of sites in the unit cell (uc) and in the 
reduced unit cell (rue) , and symmetry operations in the space 
group not including primitive translations. 



Polymer 


sites in uc 


sites in rue 


symmetries 


PAA 


6 


3 


j, 0"h, C2a, Cat., Cs 


PBA 


14 


7 


h (^h, C2, Caa, Csb 


PPR 


10 




i, ah,C2a, C26, C2c,CTvl, v2 



TABLE 11. Distortions considered for the PAA strip. For 
B distortions we identify subcases: Bi for Ai > and A2 — 0; 
and B2 for Ai = and A2 > 0. For C distortions we identify 
subcases: Ci for Ai > and Ao — 0; C2 for Ai = and 
Ao > 0; and C3 for Ai > and Ao < 0. 



TABLE IV. Distortions considered for the PPR strip, 
where j = 1,2,3 and i = 1,2,3,4. For C distortions we 
identify subcases: Ci for Ai — A4 = 0, A2 > and A3 > 0; 
and C2 for Ai = A4 = 0, A2 < and A3 > 0. For D distor- 
tions we identify subcase Di for Ai = A2 — 0, A3 > and 
A4 > 0. 



Distortion 


C2 


(Jv 


Restrictions on A; 


A 


+1 


-1 


Ai = A. = -A,, = -Ar, 
A4 = A4/ = 


B 


-1 


+1 


A, = A,^ = ~Ay = A ., 
A4 = A4/ = 


C 


-1 


-1 


A, = -Aj = A,, = -A^, 


D 


+1 


+1 


A, = Aj = A,;, = Aj, 



Distortion 


Cs 


(Jv 


Restrictions on A; 


A 


+1 


-1 


Ao = Ag = A2 = A2 = 

Ai = Aj = -Ar = -Ai, 


B 


-1 


+1 


Ao = Ag = 

Ai = -Ai = -Ai- = Ap 
A2 = -A2 


C 


-1 


-1 


Ai = -Ai = Ai/ = -Ap 
A2 = A2 = 
Ao = -Ao 



TABLE V. Ground-state Heisenberg energy per site in J 
units for the family of 7r-network polymers studied. PA stands 
for polyacetylene. The first row corresponds to the energy ob- 
tained with a single Kekule structure | K). \ vfi) stands for 
the IBR-RVB ansatz of Eqn. Q, | *3> is the 3BR-RVB 
ansatz of Eqn. (|l^, | "l>jv) is the Neel state, and | ^'jv) the 
Neel-state-based ansatz of Eqn. (|lo|). The last row corre- 
sponds to the exact ground-state energy which is known only 
for the ID case. 



TABLE III. Distortions considered for the PBA strip. All 
possible Ai, i = 1, 2, 3, 4 are assumed to be mutually indepen- 
dent. For A distortions we identify subcases: Ai for Ai > 
and A2 = A3 = A4 = 0; and A2 for Ai = A2 = 0, A3 > 
and A4 < 0. For C distortions we identify subcases: Ci for 
Ai > and the rest equal to zero; and C2 for Ai = A2 = 0, 
A3 > and A4 < 0. 



Distortion Cs ctu Restrictions on A; 



A +1 -1 A, = Aj = -A,, = -Aj, 



B -1 +1 A, = -Aj = -A,, = Aj, 



C -1 -1 A, = -Aj = A,, = -A-, 



E/JN 


PA 


PAA 


PBA 


PPR 


\K) 


-.37500 


-.37500 


-.37500 


-.37500 


1 *i> 


-.37500 


-.37500 


-.4339(3) 


-.4435(2) 


*3> 


-.41100 


-.4539(5) 






$iv) 


-.25000 


-.333(3) 


-.3214(3) 


-.32500 


*iv) 


-.4279(1) 


-.4941(0) 




-.4906(2) 


exact 


-.4431(5) 
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FIG. 1. Polymer systems. Fragments of: (a) Polyaceacene 
(PAA). (b) Poly(benz[m,n])anthracene (PBA). (c) Polyperyl- 
ene (PPR). The region between the vertical dashed lines de- 
fines the unit cell of PPR, while for PAA and PBA the reduced 
unit cell is instead identified. 



FIG. 2. Symmetry elements for: (a) PAA. (b) PBA. (c) 
PPR. 



FIG. 3. Representation of the different non-mixing Kekule 
phases of PAA, each one containing essentially one Kekule 
structure. 



FIG. 4. PAA analysis, (a) Unit cell and reduced unit cell, 
(b) Labels associated to bonds, (c) Symmetry elements cho- 
sen to label distortions: the screw axis Ca and the vertical 

plane cjv 



FIG. 5. PBA analysis, (a) Unit cell and reduced unit cell, 
(b) Bond labels, (c) Symmetry elements chosen to label dis- 
tortions: a screw axis Cs and a vertical plane (t„. 



FIG. 6. PPR analysis, (a) Unit ceU. (b) Bond labels, (c) 
Symmetry elements chosen to label distortions: a two-fold 
rotation axis C2 perpendicular to the molecular plane, and a 
vertical plane a^- 



FIG. 7. Energy as a function of the distortion parameter 
A in PAA (a) when the Neel-state-based ansatz \ ^ n) is 
considered, (b) when the 3BR-RVB ansatz \ $3) is considered. 
The curves correspond to the different distortions given in 
Table §: (□) Bi, (A) B2, (O) C2, (o) A, (.) Ci, (■) C3. 



FIG. 8. Energy as a function of the distortion parameter A 
in PBA when the IBR-RVB ansatz \ *i) is considered. The 
curves correspond to the different distortions given in Table 



[II 



(o) Ai, (□) A2, (A) C2, (O) Ci. 



FIG. 9. Energy as a function of the distortion parameter 
A in PPR: (a) when the Neel-state-based ansatz \ n) is 
considered; (b) when the IBR-RVB ansatz \ ^'i) is considered. 
The curves correspond to the different distortions given in 
Table (A) Ci, (□) C2, (O) -Di. 
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